General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




<IASA Contractor Report 156844 


Detailed Gravity Anomalies from GEOS-3 
Satellite Altimetry Data 

(V A''»-c^- 156PI*U) DFTAILFD GRAVITY ANOMALIES 
FF01 GFOS-3 SATELLITE ALTIMETRY DA^A Final 
! -port (fatten* Colusibus Labr. r Ohio.) 

“P P MC A03/MF A0 1 CSCL OflE 

G3/46 


!TM>- 12631 


Oncl as 
3 BO 1 P 






G S Gopalapillai and A G Mourad 



October 1978 


NASA 

National Aeronautics and 
Sp.»ce Admimstrat'on 

Wallops Flight Center 

Wai'ops Island Virginia .\1337 
AC 804 0?4 3411 





FORWARD 


This Investigation was performed by Battelle's Columbus 
Laboratories (BCD for the National Aeronautics and Space Administration 
(NASA), Wallops Flight Center (WFC), Wallops Island, Virginia, under 
Contract No. NAS6-2801. Nr. H. Ray Stanley served as the NASA Technical 
Monitor. 

The excellent cooperation of Professor Richard H. Rapp of the 
Ohio State University is gratefully acknowledged. As one of the GEOS-3 
Principal Investlrators, he provided all of the processed altimetry 
data in the calibration aria, most of them even before publication. He 
also made available the terrestrial gravity data and a few computer 
subroutines. Mr. James Marsh of the NAS/ /Goddard Space Flight Center 
supplied the GEM-9 potential coefficients set. 


ii 


[ # ) I 


I f k: J JJ JU-i-.l J I-i-i. -I- i— 


EXECUTIVE SUMMARY 


The objective of this investigation was to develop methods and pro- 
cedures that could be used in establishing practical applications for the high 
resolution altimeter capability of the CEOS-3 and future SEASAT satellite 
missions. The geopotential model that has direct applications in several 
areas of geophysics, mti lne geodesy, oceanography, Earth resources, etc., 
is the free air gravity anomaly. The requirements for gravity anomalies 
in potential application atcas such as oil and gas exploration are highly 
demanding in terms of accuracy and resolution. Consequently, Battelle's 
effort was directed toward development and test of a suitable technique for 
deriving mean gravity anomalies from dense altimetry data to a resolution 
that has not been achieved before. Obviously, the achievement of such a 
resolution has to be approached on a non-global basis. 

The procedure developed and applied in this study uses a combina- 
tion of both deterministic and statistical techniques. The basic mathe- 
matical model was based on the Stokes' equation which describes the 
analytical relationship between mean gravity anomalies and geold undula- 
tions at a point; this undulation is a linear function of the altimetry 
data at chat point. The overdetermined problem resulting from the 
excessive altimetry data available was solved using Least-Squares prin- 
ciples. These principles enable the simultaneous estimation of the 
associated standard deviations reflecting the Internal consistency based 
on the accuracy estimates provided for the altimetry data as well as for 
the terrestrial anomaly data. 

Using GEOS-3 data in the calibration area, several test computa- 
tions were made of the anomalies and their accuracy estimates for dif- 
ferent combinations of: 

(1) Four a priori weighting functions for anomaly 
parameters 

( 2 ) Two anomaly parameter configurations 

t3) Three data densities and distributions. 

For "profile" type or low density data, the computed anomalies 
were sensitive to the a priori weights. For such a distribution of data, 
numerical anomaly auto- and cross-covariance functions can be used as 
weighted constraints to obtain realistic estimates for anomalies and their 
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Computed anomalies for 1* x 1* blocks were compared with terrestrial 
estimates and also with a set of anomalies computed by Rapp us In* Least - 
Squares Collocation procedures. The RMS differences were 8.7 mgal and 5. 4 
mgal, respectively, which were consistent with accuracy estimates associ- 
ated with these sets (10 mgal for terrestrial data and 7.3 mgal for Rapp's 
data) . 

The accuracy estimates for 1* x 1* mean anomalies that can be 
obtained with the techniques used In this study are of the order of 2 mgal, 
which Is at lertst about four times better than the estimates for the 
terrestrial data and Rapp's data. Estimation of 30' x 30' anomalies can 
be accomplished with accuracies of about 5 mgal using the techniques em- 
ployed here. These accuracy figures are very encouraging. Indicating the 
potential applications of the altimeter data in a variety of areas. As a 
result of the lack of accurate ground truth data, proper verification of 
the fine results obtained by the procedures developed and used here Is 
not possible. Tills situation emphasizes the pressing need for a marine 
test range where detailed and accurate ground truth gravity data would be 
collected and made available to verlfv the results of GEOS-3 altimetry 
data and also of further Improved data expected of the future SEASAT 
satellite missions. 
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DETAILED GRAVITY ANOMALIES FROM GEOS- 3 SATELLITE ALTIMETRY DATA 


by 

C. S. Copalaplllal and A. C. M.urad 


1. INTRODUCTION 


This report covers activities performed by Battelle's Columbus 
Laboratories (BCL) for the National Aeronautics and Space Administration, 
Wallops Flight Center, under Contract No. NAS6-2801. The primary objective 
of these activities Is the development and testing of a suitable technique 
for estimating the fine structure of the Earth's gravity field from the 
high resolution altimeter capability of CEOS-3 and future SEASAT satel- 
lite missions. 

With the successful launch and operation of the Skylab and 
GEOS-3 satellite altimeter systems, the acquisition of detailed gravity 
data over the vast oceans of the Earth experienced a quantum jump in 
speed and detail. Those data are directly derived from the altimetry 
measurements In the form of geold undulations which describe the separa- 
tion of the geopotential level surface corresponding to the Mean Sea 

Level (MSL) and the reference ellipsoid adopted as the closest analytical 

( 1 )* 

surface representing the Earth. However, the geopotential model that 

has direct applications In several areas of geophysics, marine geodesv, 
oceanography, Earth resources, etc., Is the free air gravity anomaly. 

The advantage of the gravity anomaly model over the undulation model is 
Its higher sensitivity to the various (anomalous) features on or near the 
Earth's surface. Further, the direct observations of the gravity on the 
Earth's surface lend themselves to the natural representation of the 
anomalous geopotential field by free air gravity anomalies. 

The requirements for gravity anomalies in potential application 
areas such as oil and gas exploration are highly demanding in terms of 
accuracy and resolution. The unprecedented resolution and speed with 
which the altimeter data have become available have resulted in the possi- 
bility that the fine structure (short wavelength features) of the anomalous 


* 

References, denoted by superscript numbers, are Jt end of report in 
Section 7. 


geopotential could be determined. The determination o f this fine structure 
hae a practical disadvantage where the number of geopotential parameters 
to be determined is too large to handle, even with the largest computer 
available today. Obviously, the possibilities for solution to this pro- 
blem can only be achieved through u . of some form of non-global solution. 

Several methods and procedures can be found in the literature 
for determining gravity anomalies from altimetry data on a non-global 
basis. Some of these are purely deterministic and some are statistical, 

while others use a combination of these techniques. Examples of the 

( 2 ) 

statistical techniques are Least-Squares Collocation and Fourier trans- 
form solutions^'. A direct determination of point anomalies from alti- 

(4) 

meter data using the Inverse Stokes' equation Is an example of the 
deterministic approaches'*^'. 

The technique preferred over the above and employed in the ensuing 
study is a combination type using the analytical relationship (Stokes' 
equation) between mean gravity anomalies and geoid undulation which Is 
linearly related to the altimetry observation. This technique has several 
favorable features. The overdetermined problem resulting from the excessive 
altimetry data available can be solved using Least-Squares principles. Any 
a priori Information available on the gravity anomalies being computed can 
be incorporated simultaneously using weighted constraints. Further, any 
systematic errors Inherent In the altimetry data can be modeled and fil- 
tered out. 

This method was developed and some simulation studies made 
while the principal author was working on an Air Force contract. The 
details of that work can be found in Referenced). In the current study, 
essentially the same procedure, with minor modifications, is implemented 
with real altimetry data from GEOS-3. For the sake of completeness, a 
brief summary of the procedure is presented In Section 2, followed bv a 
short description of the various data sets (altimetry as well as terres- 
trial gravity and potential coefficients) used (Section 3). In Sections 
4 and 5, results are given for the l 8 x 1° and 30' x 30' gravity anomaly 
determinations. Conclusions from these results and recommendations toward 
any future work in this area are presented in Section 6. 
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MATHEMATICAL MODEL AMP PKOCEPUKE 

An Indicated earlier, the detail* of the technique and procedure 
used In this etudv are fully documented' 1 ' . However, for the sake of com- 
pleteness, the portion* relevant to the mathematical model and the procedure* 
for Implementing the same will he summarized here. 

J.l Mathematical Model 


For practical Implementation, the Stokes' formula, generalized to 
anv geocentric reference ellipsoid, can he written In the form: ^ 


N - N ♦ _ 

p o alii 


jEE A, su..) AO 


IJ 


u • 


( 1 ) 


where 


N - undulation at point p 

P 

N - scale constant 
o 

R - average Karth radius 
C - average gravity on Karth 
Ag - free-alr gravity anomaly 

Ao - area of the block over which Ag Is assumed constant 
4^ - distance between P and Ao 

1J - subscript identifying the block on a grldded format, 
and Is the Stokes' function given as''' : 


with 


S( V 


hx+l-5cos v^-lcos , .nix+x‘'l 




x • sin 1 




Equation ll) Implies that Ag and are compatible with respect 
to the flattening of the referet.ee ellipsoid. It also Implies that Ag be 
known In every block 1J on the Karth. However, the behavior of the kernel 
$lvl In Equation ill Is such that It suppresses the effect of the gravitv 
anomalies In the zones remote from the computation point so that only the 
lower harmonics of the gravity field In these areas contribute slgnlf leant lv 


to N 


Hie lower degree harmonics are reasonably well known and mav be 


used to advantage to compute the contribution from the zones remote from N 


In view of the Above, Che undulation N can be assumed Co be 

P 

:ompoeed of four components, N , N , N, end N. s 1 h that: 

O 1 i J 


N 

P 


N + N, 
o 1 


N, ♦ 


( 3 ) 


rhere N^, the lower harmonic contribution over the whole surface of the 
iarth, is given by: 


max n 

N .■ R T Y* (C cos m \+S sin m \) P (sin*) , 
1 , iLu nra nm no 

n«J m»0 


(4) 


ith 


(*,\) • geocentric coordinates of P 

C , S • fullv normalized potential 
nm nm 

reference field to which N 

P 


coefficients referred to the 
in Equation (l) refers 


same 


N - maximum degree to which potential coefficients are used 
nui x 

P lsin*)« fullv normalized associated Legendre function, 
ntn 

Then, N,, tb • higher harmonic contribution from a cap centered at 
he computation point is expressed as: 


N, 


4?C U Ug-Ag # )S<*)d0 


15) 


ith 


(7) 





cos m\ do 


+ S 

nm 



sin mV do , 


(b) 


here the smoothing operator 6 is set to one. The integrals in Equation u>V 

ndicate that ig is the mean gravity anomaly over the block of area A. o 
s c 

n Equation (5) implies that the integration is carried out w' spherical 
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cap of radlua Than, the higher harmonic contribution from the area out- 

side the cap Is: 


‘• c U l 


(Ag-Ag )S(<>)do 


(7) 


with o-o being the area of the remote tones, 
c 

N, would be the error In N If the actual anomaly Information In 
3 P 

the remote zones la replaced by anomalies Implied by a spherical harmonic 
potential coefficient set truncated at degree Neglecting Equation 

(3) with Equation (3) can be written In the form: 


N 

P 


+ N J J «g 

1 4i»G J a J 8 


S(tp) Ao 


( 8 ) 


where 6g ■ Ag - Ag^ , (8A) 

which can be Interpreted as the mean anomalies referred to a reference 
spherop of degree and order N 

max 

At this point we have the option of modeling any systematic errors 

that may be inherent in the altimetry observations. However, the data we 

used In this study have been preprocessed and corrected for any systematic 

( 8 ) 

errors due to tides, orbit uncertainties, etc. . Consequently, it is 

assumed that the data are free of systematic errors, other than the global 

constant parameter N , which can absorb anv bias, constant over the area 

o 

under consideration. 


2.2 Procedure For Anomaly Recovery 


The mathematical model expressed by Equation (8) forms the basin 
for the computation of mean gravity anomalies from altimeter data which, at 
thi3 point, are assumed to be in the form of geoid undulations averaged 
over the radar footprints. Every altimeter observation will result in an 
equation of rhe type (8). Conceptually, a set of such equations will give 
a meaningful solution when the number of observations equals or exceeds 
the number of unknowns (anomalies and N ) . 


t | | | I | I I I i II k r ‘ ‘ 


0 

In order to Incorporate the a priori Information, If anv, on tome 
or all unknowns, let all the unknowns also be treated as observables with 
a priori weights. Zero weights can be assigned to those unknowns for which 
no Information Is available. In this case, Equation (8) can be Interpreted 
as a condition among observables. For consideration of redundancy, the 
number of observations would be the sum of the altimetry observations and 
the unknowns with non-zero weights. The resulting redundancy In the system 
can be exploited In some type of adjustment technique. 

The mathematical model as given by Equation (8) can be rewritten 
in a more general form as follows: 

F(N, 5g, N ) - 0 , (9) 

o 

where F Is a function of the altimetry observations N, the residual anomalies 

referred to a spherop of degree N max i and the constant bias as expressed by 

Equation t8). In matrix notation, let X^ denote the vector of a priori 

estimates of the unknowns (5g, N ) and N„_ be the vector of observed values 

o Ud 

of observables (N). If the adjusted values of these quantities are, respec- 
tively, X a and N a , then: 

( 10 ) 

where V and V„. are vectors of residuals of the unknowns and observations, 
x N 

respectively. 

Linearizing Equa Ion (9), the resulting condition equation, in 
matrix notation. Is: 

V„ + AV + W - 0 , (11) 

N x 

where A is the design matrix of partial derivatives of F with respect to 

the unknowns. W is the misclosure vector obtained by evaluating F with 

the values X and N._. If the weight matrices for the observations and 
O OB 

unknowns (also considered observab les) a r e P and P , respectively, then the 

x 

Least-Squares solution for the residual vectors is: 


X - X + V 
o x 


N - N + V 
OB N 


i 


I 


j : 


7 


V - - u t pa ♦ p ) -1 a t pw 

X X 


( 12 ) 


V' - - (AV ♦ W) , (13) 

rl X 

and th« covariance matrix can be estimated from: 

I - 0 2 (A T PA + P ) _l , (14) 

X O X 

•/here o 2 la the variance of unit weight eatlmated from 
o 

0* " < v I PV u * v I p V )/» F . (15) 

O N N XXX 

with DF being the number of degrees of freedom. 

The adjusted values of the unknowns and the observations are 
obtained by substituting the vectors and In Equation (10). The ad- 
justed values of the anomalies are referred to the spherop of degree N 

max 

If the anomalies need to be referred to an ellipsoidal reference field, 
thev an be obtained by simply adding ig (Equation (b)] back [Equation 
(8A)j. 


2.3 Son-Global Solution 


The discussion thus far has been on the general feasibility of 
tecoverlng gravity anomalies from altimetry data. Hie application of the 
above method for small regions of the Earth will be discussed here. In 
fact, the Insignificant contribution to the geold undulation at any point 
from anomalies In the remote zones Is the kev to this non-global solution. 

Suppose that the anomalies within a given area, for example a 

circle of radius tl(I) centered at a point, <}, are to be determined. Then , 

the altimetrv data within an area of a concentric circle (II) of radius 

(»+-, ) are necessary to determine these anomalies so that the contribution 

Un 
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from altimetry data outside this circle la less than a predetermined value 
j . Figure 1 Illustrates this system. Further, the mathematical model 

[Equation (8)] requires that the area of anomalies be extended to a further 
distance of <|/ QS from Circle II to Circle III. This meant, there Is no 

observation equation for any observation within Circle II that will have a 
cap extending beyond Circle III. At chis point, ll should be emphasized 
that altnough all the anomalies within Circle III will be treated as un- 
knowns In this procedure, only for those within Circle I have all contributing 
factors been taken Into account. This means that the determination of only 
those anomalies In Circle I can be totally reliable. The rest of the 
anomalies will be poorly determined and their quality of determination 
will deteriorate as the distance of these anomalies from Circle I Increases. 



FIGURE 1 . SYSTEM OF OBSERVATION AND UNKNOWNS FOR 
A N0N-GL08AL SOLUTION 
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Before proceeding further, we should define the angles and if^,, 

which determine the sire of the circles. Is the sire of the cap where 

the geo Id undulations outside this cap have lnsl|tnlf leant contribution to the 

anomaly at the center of the cap to a Riven accuracy level. This cap sire 

can he determined bv considering a graph of the remote-rones contributions 

to the anomaly through the Inverse -Stoke a' formula for different cap 

sizes.' 1 ' Then, for a given error of commission due to neglect of the 

remote-tone effects, the corresponding cap sire can be Interpolated from 

the graph. Incidentally, the subscript M indicates the cap sire for 

Molodenskll's formula. Similarly, the subscript S In if^,. denotes the cap 

sire for Stokes' formula, and can also be determined In the same fashion 
(1) 
w 

2.4 Eva luatlon of Stokes' Formula 


In computing the elements of the design matrix A [Equation 111)), 
an average value of the Stokes' function as given bv Equation (2) has to be 
evaluated for every anomaly block Involved In the adjustment. Theoretical lv , 
the average value S(\f) Is expressed by: 


SU-) 


* h {l 


) do 


lib) 


Hooever, In practice. Equation lib) Is evaluated bv computing one or several 


Equation lib): 
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10* 
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ent 
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accurate 
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However, thin system of subdivision ha •* a weakness when used with altimetry 
dal*. That la, whan the computet ion point coincides with the center of the 
sub-block (l.e., p«0) , SU>> will be singular. To overcome the singularity 
in the senaltlve rattle of y, an alternative technique developed bv l.eval lols^ *' 
la used . Here, the geodetic coordinates are transformed Into Car tee lan 
coordinates using Mercator projection and the Integration required for 
averaging Sly) Is carried out analytically. The details of this technique 
can he found In Keterence (1). 

The number of times the Stokes' function Is evaluated mav be as 
manv as tens of thousands, considering the large number of altimeter 
observat Ions, the anomaly blocks Involved and the subdivision schemes used. 

It has been suggested^ 1 *' that It Is economical, with inspect to computer 
time, to do a linear Interpolation In a table of about l NO values of Sly) 
rather than to evaluate them from Equation U')> However, the Stokes' 
function Is too sensitive for such Interpo lat Ion In the range of yll)* to 20") 
Involved In the study. Tills sensitivity problem could be overcome bv 
modi tying the procedure as follows. 

The function Fly) defined as: 


Flp) 


sin 


Sip) 


117) 


(o) 

Is known to be stable for values of y within the range undei const do rat Ion. 
Therefore, the values of Ely), Instead of Sly), are tabulated. Then the 
value of Fu l for a given value of y Is obtained bv linearly Interpolating 
between the tabulated values, and the corresponding Sly) Is slmplv obtained 
using Equation l l 7). 

The values of Fly) are tabulated as follows, with the values of 
in degrees: 

Foi 0,1 » i* ' 1 at everv hundreth of a degree 
1 < p < 5 at every tenth of a degree 
y ■* S .it every degree 

ill) 

y < 0.1 a constant value of 1.000’ Is used 
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3. DATA DESCRIPTION 


The various types of data used in this study are described here. 
These include: 

(1) Altimetry data 

(2) Terrestrial anomaly data 

(3) Spherical harmonic potential coefficients set 

(4) Numerical auto- and cross-covariances between 

1* x 1* and 5* x 5* anomaly blocks with distance <(». 

Before proceeding further with a description of each data set, it is neces- 
sary to define the area selected for the anoiaaly recovery. The basic cri- 
terion for the location of this area is the coverage of the altimetry data 
available at the time the computer programs were being developed and tested. 

This preprocessed data, free of systematic bias and orbit errors to within 

( 12 ) 

78 cm, were made available to us by Professor Rapp of The Ohio State 

University. The best coverage of these data was in a 5® x 5® equal-area 

block whose latitude and longitude bounds are 40®-35® N and 291®-297® E, 

respectively. This block, which is approximately equivalent to a cap of 

radius iji - 2?5 (Figure 1), is within the calibration area. 

For in Figure 1, a value of 3® is chosen, which corresponds 
UM 

to a truncation error of about 1.5 mgals on the anomalies to be computed. 
This value is taken from the graph on page 23 of Reference (1) and cor- 
responds to Standard Earth (SE)II^^ potential coefficients truncated at 

(14) 

degree 16. Considering the improved Goddard Earth Model iGEM) 9 used 

in this study, this error of 1.5 mgals should be a very conservative 
estimate. 

In the case of the truncation cap angle, q> 0 » for the Stokes' 
formula, a value of 22 degrees which corresponds to a truncation error of 
Just over 0.6 m is chosen.^ Once again, this error corresponds to the 
SE II potential coefficients truncated at degree 16. Consequently, the 
outer circle (III) is of radius 27®5; the radius of the middle circle is 
5 ? 5 a; ^ that of the inner circle is 2?5. Assumirg for practical purposes 
square areas Instead of circular caps, terrestriil anomalies are required 
over a 55® x 55° block and the altimetry data would be required over a 
11® x 11® block. 


I t I I I V I 1 II -i ill* ( .» > 


12 


3.1 Altlme try Pat a 


These data are taken from a set of sea surface topography heights 

( 8 ) 

ard their standard errors in the CEOS-3 calibration area, supplied by Rapp' 
The data are derived from the raw altimeter data corrected for tides and 

( 8 ) 

are free of bias and orbit errors to about 65 cm In the calibration area ' . 
There are about 15,000 observations within the 11* x 11* area selected for 
this study. However, the maximum number of observations used in this study 
is 7475, taking every other observation. The geographical distribution 
of these data is as shown in Figure 2. This data set also had geoid un- 
dulations implied by the GEM-9 potential coefficients set truncated at 
degree 20. 


3. 2 Terrestrial Cravlty Data 


The anomaly data used in this study consist of two sets. The 
5° x 5° equal-area anomalies and their standard errors were taken from 
Appendix A of Reference (7); the 1° x 1° (approximately) equal-area anomalies 
and their standard errors were supplied by Rapp from his tape called 
"August 1976 Tape"^ 7 \ These anomalies are referred to the Geodetic Reference 
System 1967^^. 


3.3 Spherical Harmonic Potential Coefficients 


These data are a set of fully normalized potential coefficients 

( 14) 

complete to degree and order 20 from the GEM-9 and referred to the 
Geodetic Reference System 1967^^. 


3. 4 Numerical Autocovariance Model 


The sizes of the blocks of anomalies involved in the recovery of 
1° x 1° anomalies are 5 and 1 degrees. The autocovariance models for 1° x 1° 
anomalies and 5° x 5° anomalies and their cross-covar lances are computed 
using the subroutine CGVA^^\ which uses the assumption that bv varving 
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the heights of the points P end Q In the point anomaly autocovariance 
function C(P,Q) points and Q., can be found, for which anomaly covariance 
C(Q^,Q,) gives a good approximation to any (sav 1* x 1*) mean anomaly 
covariance function. For example, • Q, - 10.4 km will give covariances 
for 1* x 1* mean anomalies. Similarly, 5* x 5* mean anomaly covariances 
can be approximated by C(Q^, for • 98.45 km l * h \ The cross-cover lance 

between the 1* x 1* and the 5* x 5* mean anomalies are computed bv setting 
• 10.4 km and Q, • 98.45 km. A set of these covariances and cross- 
covariances computed by COVA for various angular distances, are given 
In Table 1 for Nmax • 20. 


TABLE 1. MEAN ANOMALY COVARIANCES AND CROSS- 
COVARIANCES FOR 1-DEC AND 5 -DEC BLOCKS 




lxl Dee 

1x5 Deg 
‘J5AAllL_ 

5x5 Deg 
^ragal-) 


0.00 

654.41 

208.81 

103.64 

.50 

433.78 

190.95 

98.83 

1.00 

268.48 

152.46 

86.28 

1.50 

170.20 

112.55 

69.89 

2.00 

105.90 

78.30 

32.95 

2.50 

61.43 

50.72 

37.31 

3.00 

29.83 

29.10 

23.78 

3.50 

7.17 

12.47 

12.57 

4.00 

-8.95 

-.06 

3.59 

6.00 

-33.78 

-22.18 

-14.31 

8.00 

-26.93 

-19.60 

-14.08 

10.00 

-10.69 

-8.55 

-6. 70 

12.00 

3.54 

1.90 

.90 

14.00 

10.72 

7.52 

5.24 

16.00 

10.48 

7.68 

5.61 

18.00 

5.41 

4.15 

3. 17 

20.00 

-.81 

-.38 

-.13 

22.00 

-5.09 

-3.61 

-2.55 

24.00 

-5.98 

-4.37 

-3.19 

26.00 

-3.84 

-2.89 

-2.16 

28.00 

-.29 

-.31 

-.29 

30.00 

2.74 

1.94 

1.38 

35.00 

2.02 

1.53 

1.15 

•« 0.00 

-2.76 

-2.01 

-1.46 

45.00 

-.17 

-.16 

-.14 

50.00 

2.13 

1.57 

1.15 
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£. 1 * x 1* MEAN GRAVITY ANOMALY RECOVERY 

Since this study la a continuation of that reported In Reference (1), 
a brief auonary of the previous results and conclusions may be helpful before 
presenting the results of the current study. The previous study was based 
purely on simulations where there were only one altimeter observation per 
block (l*xl*) over a 30* x 30* area. Mean anomalies over 1* x 1* blocks 
were recovered from these observations In a 10* x 10* area. 

The highlights of the results were that: 

(1) The recovery of 1* x 1* mean gravity anomalies 
was feasible with the procedure outlined earlier 
In this report 

(2) A posteriori estimates for standard deviations 
obtained from the density of observations used 
(one per block) were about 19 mgal. 

(3) Exclusion of the anomalies in the outer zone 
(Figure 1) as unknowns resulted In modeling 
errors (aliasing effect) being Introduced In 

the anomalies being recovered; however, the aliasing 
effect on the a posteriori standard deviations was 
not very significant 

(4) When these anomalies were Included, the system of 
equations became unstable so that (a priori) wieghted con- 
straints were necessary for anv realistic results 

(5) Minimum requirement for the density of observations 
was one per block. 

These results form the basis for this study. 

The major difference between the two studies Is that the current 
study uses very dense real data while the previous study was based on simu- 
lation data which were very sparse. Consequently, the present study will 
attempt to examine several Issues raised at the conclusion of the last study. 
The Principal Issues Include: 

(1) The performance of the technique and procedure 
with real data 

The effect of the high density data on accuracy 
estimates of the anomalies recovered 


( 2 ) 
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(3) The use of swan anomalies of larger blocks In the 
outer and middle cones (Figure 1) 

(4) The use of numerical auto- and cross-covariances 
for weighted constraints on the anomaly parameters. 

It was hoped that answers to these Issues would help In the formulation 
of an optimum procedure for improving our knowledge of the geopotentlal 
field using satellite altimeter data. 

4.1 Mean Anomaly Parameter Systems 

With the cap sizes chosen (in Section 3) for a non-global solution, 
the parameters will Include all the equal-area mean anomalies within the 
area of a 55* x 55* equal-area block. This will result In 3025 1* x 1* 
equal-area mean anomalies, which is too large and expensive to handle for 
several tests to be performed In this study. Consequently, 5* x 5* equal- 
area mean anomalies are considered for the outermost 20* of the outer zone 
and 1* x 1* anomalies for the rest of the area, as shown in Figure 3, where 
the extent of the altimeter data Is shown with the rectangle with broken 
lines. As a result, the total number cf anomaly parameters would be 337 
(112 - 5* and 225 - 1°). In the other tests, 5° anomalies are considered 
for both the middle and outer zones, In which case the number of anomaly 
parameters Is 145 (120 - 5® and 25 - 1°). In the ensuing discussions, the 
system with 337 anomaly parameters will be referred to as System A, and 
System B would be that with 145 anomalies. It should be noted that there 
would be one error model parameter, N , to be determined In addition to 
the anomaly parameters. 

4.2 A Priori Weights For The Parameters 

The a priori weight matrix In Equation (12) can be considered 
to consist of two submatrices, given by: 



FIGURE 3. SYSTEM OF EQUAL AREA MEAN ANOMALY PARAMETERS 
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(18) 


Implying Chat correlation between the realdual anomallea fig and N Q are zero. 

There are eeveral waya the weight matrix for the nnonaly 
parameters can be defined, depending on che type of a priori Information 
available on than. In the ensuing test for the mean anomaly connotations, 
Pggls defined In four ways: 

(1) It la set to zero (P. • 0), Implying that we have 

6g 

no Information on them. 

(2) Using the estimates provided by Rspp^ for the standard 
errors of the anomaly values determined from terrestrial 
measurements, Pj R ls defined as a diagonal matrix D 1 
where the elements of D are the square of the standard 
errors. That Is: 



(19) 


implying that the anomaly estimates are Independent of 
each other. 

(3) Assuming that the anomalies are stochastic quantities, 
their weight matrix is defined using the auto- and 
cross-covariances among them Interpolated from the 
values presented in Table 1. If the auto- and cross- 
covariance matrix is C, then: 



( 20 ) 


(4) If the terrestrial estimates are assumed to consist of 
signals and noise^ 7 \ then: 


P 5g - (C+D)" 1 , (21) 

where C is the auto- and cross-covariance matrix and 
D is a error covariance matrix assumed diagonal as 
in Item (2) above. 
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For the purpose of determining the e priori weight for the 
constant, N 0 , term, any error In la considered to come from two major 
sources: 

(a) Uncertainty In the aemlmajor axle of the mean Earth 
ellipsoid 

(b) Uncertainty in the equatorial gravity. 

(18) 

The accuracy astlaata for the aemlmajor axle le given by Rapp as 

2.5 m. The error contribution from an uncertainty of 1 mgal In the 

(9) 

equations! gravity for a cap of 20* la about 2.6 m. Therefore, the error 
due to both sources la 3.6 m, which, rounded to 4 o la u.ied for the accuracy 
of N c -term. Consequently, 

P N - 1/16 m‘ 2 . (22) 

o 

For a priori values of anomalies, the terrestilal data described 
In Section 3.2 are used when P^ la defined according to Equations (19) and 

(21). Satellite anomalies as defined by Equation (6) are used when P. Is 

-1 * 8 
either 0 or C . However, lr. all the ensuing computations, the a priori 

value for N la assumed to be zero, 
o 

4.3 Tests With "Profile" Altimetrv Data 


At the time of the computer programs being modified/extended to 

accommodate the real data, it was decided to test them with the only pre- 

( 12 ) 

processed altimetry data available in the literature . These data, which 
will be referred to as the profile data, were supplied to us on computer cards 
by Rapp. The distribution of these data (In the calibration area) is shown 
on Figure 4. It can be seen from this figure that only five passes are 
within the block where the gravity anomalies are sought. 

The system of gravity anomalies assumed for these tests Is Svstem 
B, shown on Figure 3, where the total number of unknowns Is 146. Four tests 
were carried out with the weight matrices P, as described in the last 
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FIGURE 4. DISTRIBUTION OF THE TEST ("PROFILE") ALTTMFTFR IVTA 
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In the first test where P. was set to zero, the system of normal 

5g 

equations was unstable. The results of the other three tests are sum- 
marized in Table 2, where the blocks In which there are no altimeter data 
available are excluded. In this table, the average a priori standard 
deviation and the average a posteriori standard deviation are self-ex- 
planatory. The Root Mean Square (RMS) anomaly difference refers to the 
difference between the terrestrial anomaly estimates and those computed In 
these tustb. The results of Table 2 show that there Is very little dif- 
ference between the three weighting procedures. However, comparing the 
recovered anomalies In each of these tests, we find that while the RMS 


Is only about 


difference between the tests where P. ■ (OD)”* and P. - 

6g <g 

1 *l«l, It Is about 12 ragal for the test where P^ • D * and each of the 

other two. This Indicates the strong possibility that the a priori relative 

weights for the anomalies which are computed from terrestrial estimates 

are Incorrect. 

How much of an Improvement these anomalies obtain In their esti- 
mates from the altimeter data Is demonstrated In Figure 4, where the a 


posteriori standard deviations 


(16) 


corresponding to the case of P 


Jr 


,-i 

c 


are given. These numbers show that the anomalies of blocks where there are 
no altimeter data are very poorly determined. Their determination Is best 
when the altimeter pass goes through the center of the block or when the 
data are dense and uniformly distributed within the block. 


TABLE 2. RESULTS OF ANOM^^Y RECOVERY 
WITH "PROFILE" DATA 


Statistical Parameter 

% - (C+D)-> 

P. - c- 1 

P. - D- 1 

Average standard deviation (lo), rngal 
A priori 
A posteriori 

28.1 

9.1 

26.3 

8.4 

10.0 

7.2 

RMS difference with terrestrial 
anomaly, mgal 

16.3 

16.8 

16.4 

Variance of unit weight (o ) 

o 

0.8b 

0.85 

1 .08 

Constant N , m 

0 

-7.0 ♦ 1.4 

-5.1 t 1.2 

-2.1 0.7 


22 


4.4 Tests With Low Penalty Altimeter Data 


The tern low density refera to the set of 1496 data points which 
la one of every five data points shown on Figure 2. The difference between 
this and the profile data as discussed in the last section Is that the low- 
density data are more evenly distributed In the area than the profile set. 

In order to study the effects of these low-density data on this anomaly 
recovery, three nore sets of computations were made with the same System (B) 
of anomaly unknowns. Since the difference between using P^ ■ C 1 and • 
(OD) 1 Is small, as seen from the previous tests. It was decided that the 
three new determl nation* would use: 

<•> P A» " 0 

-1 

(b) P. - (C+0) 

4g -1 

to >„■»' . 

The results of these determinations are summarized In Table 3. 


TABLE 3. RESULTS OF ANOMALY RECOVERY 
WITH LOW DENSITY DATA 


Statistical Parameter 

P. - 0 

fig 

P. - (C+D)~ l 
fig 

P. - D" 1 
fig 

Average standard deviation (lo), mgal 
A priori 
A posteriori 

00 

5./ 

28.1 

6.0 

10.0 

6.1 

RMS difference with terrestrial 
anomaly, mgal 

14.8 

13.2 

13.0 

Variance of unit weight (o q ) 

1.06 

1.24 

1.62 

Constant N , m 
o 

0.7 ♦ 4.1 

-3.6 i 2.3 

0.5 t 1.0 


Once again, on the average, there is very little difference between 
results of toe various weighting procedures. Results in Table 3 indicate 
that the a posteriori standard deviation estimates as well as the RMS dif- 
ference between the computed and terrestrial anomalies have improved somewhat over 
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those with the profile data. Also, the RMS differences between the anomalies 
determined with P^ • (C+Dl ' and • D ' have dropped to 7.5 mital compared 
to 12 mgal with the profile data. In fact, the difference Is less than 5 
mgal for 19 of the 25 anomalies computed. However, the reasons for the 
differences of greater than 10 mgal In the other six anomalies are not 
clear at this time. It should also be noted that the svstem of normal 
equations appears to be stable for P^ * 0 for the anomalv and altlmetrv 
data systems under consideration. 


4.5 Tests With High Density Data 


In these teats, the altimeter data used are as shown on Figure 2, 
(7475 data points! which la five times denser than the low-density data 
discussed In the last subsection. These tests are aimed at examining the 
effect of hlgh-densltv data on the accuracy estimates of the anomalies being 
determined. Here again we use the same weighting functions and the anomalv 
parameter System B that were used with the low-denaltv data. The results 
corresponding to these tests are given In Table 4. 


TABLE 4. RESULTS OF ANOMALY RECOVERY 
WITH RICH DENSITY DATA 


Statistical Parameter 

P. - 0 

-sr 

P, - ( C + D )" 1 

-sr 

p - rr 1 

-'R 

Average standard deviation (lo>, mgal 
A priori 

a> 

28.1 

10.0 

A posteriori 

2.b 

2.5 

2.7 

RMS difference with terrestrial 

anomaly, mgal 

14.8 

14. \ 

14.1 

Variance of unit weight (o 1 

o 

l .Ob 

1.13 

l . 34 

Constant N , m 

2.4 ♦ 1.5 

4.8 ♦ 1.9 

0.17 • 0.8 

o 

— 
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This la Again another example of the consistent agreement among 
the v«i 1 ous weighting algot'lthma. rite dominance of the altimeter data can 
he aeen from agreement between the computed anomallea under the varloua 
weighting procedurea. Hie RMS difference for the caaea P^, ■ 0 and P^ • 
(OO) * la 8.6 mgal. For the caaea P^, - (OD) * and P^ • D \ thla dif- 
ference la only 4 mgal. Comparing the resulta of the hlgh-denaltv data with 
thoae of the low-denaltv data, the estimated standard deviation averaged for 
the 25 blocks Is Just under half fo. five times Increase In the density of 
altimeter data. The RMS difference In the computed anomalies la only 4.6 mga 

4 .6 Testa With 1* x 1* Anomalies of System A 

The anomaly System A (Figure J) consists of 1° x 1" anomalies 
extending 2 degrees beyond the middle zone (Figure 1). The total number of 
anomaly parameters in this svstem Is 33'’. of which 225 are 1" x l* anomalies. 
The number of observations corresponding to the htgh-dcnslly data Is 7475. 
Computer timewise. It would be verv expensive to use t he hlgh-densltv data 
with the svstem of 318 unknowns Including th«. N -term. Consequently, we 
decided to examine the effect of extending the 1* x l® anomalies beyond the 
middle zone with the low-density data which has 14Sn data points. The 
weighting functions remained the same. The results are summarized In Table 5 


TABLE 5. RESULTS OF ANOMALY RECOVERY 
FOR ANOMALY SYSTEM A 


Statistical Parameter 

P, - 0 

Jr 

P. - (C+0)- 

Jr 

l p - p-1 

Jr 

Average standard deviation (lo), mgal 




A priori 

>v 

28.1 

10.0 

A posteriori 

23.4 

4 • 4 


RMS difference with terrestrial 
anomaly, mgal 

ld.0 

11.7 

10.5 

Variance ot unit weight (o ) 

o 

0.72 

0.82 

0 . ') 3 

lonst.int N , m 
o 

0.1 t 2.d 

— « . 8 ♦ l . h 

-7.8 t 0.6 
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Unlike with Che anomaly System B, the system of normal equations 
Is un uable for zero a priori weights on the anomalies. The RMS difference 
between the computed and the terrestrial anomalies Is the lowest (about 
11 mgal) among all the different sets of computations we have done In this 
study. 

The effect of extending the smaller anomaly blocks (1* x 1*) 
beyond, the middle zone can be best seen In a comparison In Table b between 
the results obtained w<th Systems A and B. The weighting function used In 
this comparison Is (C*D) * . 

TABLE 6. COMPARISON OF RESULTS WITH SYSTEMS A AND B 


Statistical Parameter 

System A 

System B 

Average a posteriori standard deviation (lo), mgal 

4 . 4 

b.O 

RMS difference with terrestrial anomaly, mgal 

11.7 

13.2 

RMS difference with P • D~ l , mgal 

5.8 

7.5 

Variance of unit weight (o ) 

o 

0.82 

1.24 

Constant N , m 
o 

-4.8 t l.b 

-J.b ♦ 2.3 

RMS difference (A-B> 

9.8 

mgal 


The results In this table demonstrate the overall Improvement In 
performance with System A when compared with Svstem R. Yet, the differences 
are not dramatic. However, the anomalies recovered In both the systems are 
significantly different iRMS difference ™ d.g mgal) but still consistent with 
the accuracy estimates associated with the anomalies. Most of this dif- 
ference could be attributed to the aliasing effect resulting from approxi- 
mating the detailed anomalous gravity structure in the middle zone with 
larger block anomalies. 

The accuracy estimates (4.4 mgal) obtained with low-density data 
are excellent considering the accuracy of the terrestrial data available 


today. This estimate oar. be further Improved to about 2 - 2 . 5 rngal with the 
hlgh-d*n*lty data, as we have seen with the result* with the anomaly 
System I' (Table 4). 


s.7 Co mparison With Rapp’ s A noma 1 tea 

In order to determine how good these computed anomalies are, we 
provide the following comparison In Figure 5. Four candidate sets of 
anomalies are selected foi this comparison: 

(1) Set l - Tills set corresponds to Svstem A, as 

discussed In the last subsection, with a 

priori weights P, - (C+Pl In this svstem, 

eg 

1* x 1* anomaly unknowns extended 2 * bevond 
the middle zone. 

(2) Set 2 - Tills correspond# to System U with 5* x V 

anomalies used *or the middle and outer tones, 
and high density altimetry data. Again, the 
weighting function Is (t>P> 

l" Rapp's anomalies - computed using Least -Squares collocation 

i >> i 

techniques . These anomalies correspond to 

1° x l* equiangular blocks. The difference be- 
tween the equiangular and equal-area blocks In 
the area under conslderat Ion Is that two equi- 
angular blocks term an equal-area block within 
a 2 * longitudinal band In the middle of the 
*>* x *>* equal-area block, as shown In Figure S. 

In such cases, equal-area anomalies are derived 
from averaging the equiangular anomalies and 
accuracy estimates. 

(■*1 Terrestrial anomalies - as described In Subsection 1.2. 
Quantitative dltterences among the various at.om.ilv sets compared 
are computed in the form of RMS differences and shown In Table " . The larger 
differences occur with Set J and the smallest one Is between Rapp's and 
ter estrlal anomalies followed closely bv Set 1 Rapp differences. In a 
further look at the differences between Set 1 and Rapp's anomalies. It was 
found that significantly large differences occur In the blocks where the 
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FIGURE 5. COMPARISON OF COMPUTED ANOMALIES 
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equiangular anomalies had to be averaged to get equal-area anomalies. Con- 
sequentlv, we compared the two sets again excluding these five blocks, and 
the agreement Is significantly improved (5.4 mgal). Similar comparison 
between Set 1 and terrestrial anomalies shows the aame Improvement (8.7 mgal). 

The larger differences in these five blocks mav perhaps bo due to the fact 
that the equal-*area blocks are not equal area at all and the equiangular 
blocks are, In fact, equal area In a nou-global solution. Consequently, 
the l* x 1* covariance function Is Improperly used for 1* x 2* blocks. How- 
ever, this effect Is considerably reduced by the dominance of the hlgh- 
denslty altimeter data as we have seen with System H. 

An accuracy estimate of 2.5 mgal was obtained for 1® x 1* equal- 

area anomalies with the hlgh-densltv altiraetrv data. The corresponding 

( 19) 

estimates for terrestrial data Is 10 mgal and 7.3 mgal for Rapp's 


TABl.E 7. RMS DIFFERENCES BETWEEN THE ANOMALY SETS 


RMS Differences 


Anomaly Set (mgal) 

Set l - Set 2 q.g 
Set 1 - Rapp g , 7 
Set l - terrestrial 11.7 
Set 2 - Rapp 11.4 
Set 2 - terrestrial 13.2 
Rapp - terrestrial 7.5 
Rapp - Set l (20 equiangular blocks only) 5.4 
Set 1 - terrestrial (20 equiangular blocks only) 8.7 
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4.8 S umme-y of 1* x 1** Anomaly Recovery 

As outlined at the outset of this chapter, ve examined three 
aspects of the determination of gravity anomalies from real altimetry 
data: 

(1) The effect of data distribution on the accuracy of 
the anomalies determined 

(2) The effect of substituting smaller blocks for the 
larger blocks within the middle and outer zones 

(3) Use of auto- and cross-covariances for gravity 
anomalies as a priori weighting functions. 

The results can be summarized as follows: 

(1) It ha. been demonstrated that the technique and 
procedures used here can be used effectively with 
real altimetry data in the recovery of 1° x 1* mean 
gravity anomalies. 

(2) Realistic estimates for these anomalies can be 
determined with varying degrees of accuracy depending 
on the density and distribution of the altimetry data 

(3) When the density of observation is low or if only a 
few profiles of altimeter data are available, the use 
of auto- and cross-covariance functions as a priori 
weights can result in improved anomaly estimates. 
However, the anomaly blocks should have at least one 
profile passing through for any significant improvement 
in the anomaly estimates from altimetry data. 

(4) The use of larger anomaly blocks in the middle and 
outer zones to reduce the number of parameters 
results in significantly large aliasing effect/model 
errors introduced in the anomalies being sought within 
the inner zone. However, the use of 1° x 1° blocks 

to at least 2 to 3 degrees beyond the middle zone and 
0 ° x 5° blocks in the rest of the out r zone may be 
optimum in the sense that the number of parameters would 
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be small enough for practical purposes ami at the 
same time the resulting aliasing effect would be small 
enough to be neglected. 

(5) 1* x l' mean anomalies can be recovered to about 4.4 

mgal accuracy with low-density data (one of every five 
observations shown In Figure 2). This accuracy can be 
Increased to 2-2.3 mgal with hlgh-densltv data (all 
data points shown on Figure 2). 

(b) The anomaly set considered to be the most realistic of 
those computed In this studv compares verv well with 
Rapp's estimates and with the terrestrial estimates. 

The accuracv estimates obtained both with low- and hlgh- 
densltv data show considerable Improvement over Rapp's 
as well as those of the terrestrial estimates. 

These comparisons, however, are meaningless unless we have 
ground truth data with an accuracv better than the 2-2.5 mgal accuracy 
claimed In this study to verify the validity of these claims. Until then, 
we can only sav that the results obtained In this studv are comparable with. 
If not superior to, the best estimates available todav. 


I 
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3. 30* x 30' MEAN ANOMALY RECOVERY 


IC may be recalled Chat the objective of this study is to develop 
practical applications for the high resolution altimeter cspabllity of 
GEOS-3 and future SEASAT satellite missions, with the primary effort being 
the development and test of suitable techniques for gravity anomaly deter- 
mination in the marine area which is one of such applications. We have 
effectively demonstrated the Least-Squares technique In recovering 1* x 1* 
equal-area mean anomalies from both the low and high density altimetry data. 
Now, It remains to be seen whether this technique can be used to resolve 
finer structures of the gravity field from the high density altimeter data. 
In the ensuing discussions, we will present some preliminary results of a 
30' x 30' mean anomaly determination. 

5.1 Parameter System Fot The 30' x 30' Anomaly Recovery 

The high density data were used in determining the 30' x 30' 
mean anomalies. In order to keep the number of data points and the number 
of parameters in the system low for economical reasons, we decided to re- 
cover 36 30' x 30' approximately equal-area mean anomalies In a 3° x 4® 
block bounded by latitudes 36®-39® north and longitudes 292*-296® east 
(Figure 6). These blocks were subdivided from nine 1“ x 1° equal-area 
(approximate) blocks which were used in the 1® x 1® anomaly recovery. High 
density altimeter data from a 5* x 5® equal-area block (1858 data points) 
are used for this computation. As can be seen from Figure 6, these data 
extend 1 degree beyond the inner zone in which the anomalies are sought. 

In the middle zone, which is a band of 1 degree around the inner 
zone, 1® x 1® anomalies are used with another 20-degree band of 5® x 5® 
anomalies In the outer zone. This system results in 132 anomaly parameters 
and one error model parameter N 0 . 

5.2 A Priori Information On The Anomalies 


Since terrestrial estimates for the 30' x 30' blocks were not 
readily available, the terrestrial value of the 1® x 1® block of which the 
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30' x 30' blocks were subdivisions were used for sll Che four subdivisions. 
Tsrrsstrlsl estimates were used for chs 1* x 1* snd 5* x 3* blocks. 

Equation (19), P tf ■ D* 1 «od«l, was usad for weighting the anom- 
alies where the elements corresponding to the 30' x 30' blocks were set to 
zero. Auto- and crosa-covarlance models were not used since they were not 

readily available. However, an approximate model could be computed with 

(16) 


the subroutine COVA 
(16) 


Rapp 


employing the approach used by Tschernlng and 
Th'.s approach uses the similarities between the smoothing 


operator and 


« Q / 


1+2 


, where R^ Is the average Earth radius and h^, 


the height of a point above the Earth such that a point covariance function 
C(Q,Q) will adequately approximate the covariance function of x 6* block 

anomalies . 

B, can be expressed In terms of the radius of a circular cap 

as: 


1-cos 6 


— \v 

21+1 1 - 


.(coa ) - P , . ( cos * ) 
1 o t+l o 


(23) 


where P. Is the Legendre Polynomial of degree i. Now, can be approxi- 

* o 

rt° 

mated equating the area of circular and square caps; that is, . 

o | w 

For 30' x 30' blocks, 6 ■ 0 ? 5 . Thus, knowing the values of 6., a value of 
h^ could be found such that: 


ii+: 


R +h n 

e Qi 


(24) 


for all 1. Then h^ can be used as an input for the subroutine COVA to 
compute the numerical auto- and cross-covariance model. Because of 
limitations of time and resources In this project, however, computations 
of these covariances could not be accomplished. 
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5.3 Results and Discussion 

The 30' x 30* mean anomalies as computed with high danalty aid- 
Htir data art prasentad In Flgura 6. Their avarage standard deviation 
(Id) was 5 mgal with the estimate of the variance of unit weight being 0.6. 

The error model parameter N waa estimated to be - 1.76 t 0.5 meters. 

o 

These numbers look very realistic except In the upper right 
corner of Figure 6, where unusually high and low(-69 and 28) values arc 
noted In the adjoining blocks. A closer look at Figure 2 Indicates that 
this block has poor "across the track" data coverage. These values mav be 
due to Chat; on the other hand, they mav be reasonable. Only good ground 
truth data can resolve such questions. 

As we mentioned earlier, the terrestrial estimates for these 
anomalies were not readily available for comparison at the time this 
report was written. However, 1* x 1* mean anomalies are averaged from 
these 30' x 30' anomalies and compared with those from Figure 5 correspond- 
ing to Set 1, Rapp's, and the terrestrial sets. The results of this com- 
parison Is provided In Table 8 in the form of RMS differences. In the 
comparison with Set 1, the larger blocks, where unusually larger differences 
due to Improper covariance function were noted, were excluded. 


TABLE 8. COMPARISON WITH OTHER ANOMALY SETS 


Averaged 1“ x 1* Anomalv MinuB 

RMS Differences 
^mgal) 

Terrestrial 

6.7 

Rapp 

4 .4 

Set 1 (Figure 5) 

4.8 


These comparisons are very good considering the standard devi- 


ations associated with them. However, this agreement is not an indication 


I 


that the Individual 30' x 30* mean anomalies obtained In these computet Ion* 
are resolved with good accuracy, but It certain'.? Indicate* that the esti- 
mate* arc very realistic. Accurate ground truth data are required to 
verify these estimates. 

The accuracy of 5 agal for the 30' x 30' anomalies Is excellent 
considering the data distribution, which Is poor In the direction across 
the track (satellite pass). It would be Informative to see how the more 
uniformly distributed data expected from the SEASAT missions will affect 
these accuracv estimates. 


6. CONCLUSIONS AND RECOMMENDATIONS 


The results of the study presented in Sections 5 and 6 of this 
report led to the following conclusions: 

(1) The technique and procedures developed and used 

in this study are very effective in the determination 
of detailed mean gravity anomalies to the extent of 
1* x 1* and 30' x 30' blocks. 

(2) The accuracy estimates for the computed anomalies are 
highly dependent on the density and distribution of 
data points. 

(3) The anomalies computed are sensitive to a priori 
relative weights us«:d as constraints. This sensitivity 
reduces with Increased density and uniform distrl- 
butlon of the data points. 

(<•) The use of the accuracy estimates given for the 

terrestrial anomalies for a priori weight constraints 
results In poor relative weights among the anomalies. 

(5) The use of auto- and cross-co .ariance models for 
anomalies as a priori weight constraints Is very 
effective In determining realistic anomaly estimates 
when the density of data points is low or if only a 

few profiles of altimetry data are available In the area. 

(6) A block should have at least one profile passing 
through it for any significant improvement in the 
estimate of its value. 

(7) Approximation of small blocks by larger blocks in the 
middle and outer zones reduces the unknowns and in- 
creases the stability of the normal equation system. 
However, this introduces significant model errors/ 
"aliasing" effect on the anomalies to be determined. As 
a compromise, the use of the small blocks in at least 
the middle zone results in a reasonable size of the 
parameter system and small modeling errors. In fact, 
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the small blocks were extended 2 degrees beyond the 
middle and Into the outer rone In this study. 

(S) When the auto- and cross -covariance* are used for a 
priori weighted constraints. Ideally the blocks 
where the anomalies are determined have to be truly 
equal In area. In a non-global solution, equiangular 
blocks would be approximately equal area. The equal- 
area subdivision, as used In this study, where some 
blocks are double the area of the others, results In 
Improper relative weights. However, this needs to be 
examined further. 

(9) With the low density data (149h data points over a 
11* x 12* block), 1” x 1* mean anomalies can be de- 
termined with an accuracy of about 4 mgal. Indications 
are that this figure can be Improved to about 2 mgal 
with high density data (7475 data points). A prelim- 
inary determination showed that 30' x 30' mean 
anomalies can be estimated with an accuracy of about 
5 mgal with the high density data. 

(10) Computed anomalies compared well with those computed by 
Rapp using Least-Squares Collocation techniques and with 
terrestrial estimates; the accuracy estimates obtained 
In this study are slgnlf leant lv better. However, due 
to the lack of accurate ground truth data, we cannot verlfv 
how good each of these determinations is. 

Even though significant results have been obtained considering 
the available time and resources In this studv, there are several areas 
tnat need further examination. Some of these are listed below: 

(1) Some of the equal-area blocks used here are larger 
than the other blocks due to the subdivision scheme 
used. The 1® x 1° autocovariances used for a priori 
weight constraints appear to have unfavorable effects on 
the anomaly estimates. These effects need to be examined 
by using either truly equal-area blocks or equiangular 
h locks . 


I » * • • • « • 
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(2) The 1* x 1* anomalies must be recomputed using high 
density data with Parameter System A, where the 1* x 1* 
anomaly parameters are extended 2 s beyond the middle 
tone. This computation was not done In this utudy due 
to the limitations of time and resources. 

(3) The data distribution In Figure 2 shows poor cross- 
track coverage. It would be Instructive to see the 
affect of more uniform data coverage on the computed 
anomalies and their accuracy estimates. 

(4) The anomaly parameter svstem used for the 30' x 30' 
anomaly determinations must be reexamined to see 
whether anv improvement to the anomaly and their 
accuracy estimates could be made. 

(5) The possibility of resolving mean anomalies over 

still smaller blocks (e.g., 15' x 15') must be examined, 

(6) However Improved the anomaly and the accuracy 
estimates mav be, they are meaningless unless we can 
verify them with proper ground truth data. The ground 
truth data presently available are totally inadequate 
with respect to coverage and accuracy. There should 

be a concerted effort on tne part of the agencies respon- 
sible for geodesy to establish a marine ca 1 lb -at ion 
range where proper ground truth data could be established 
for verifying the high resolution altimeter data and 
their derivatives. 
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